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Abstract 

The collective flow of nucleons and that of fragments in the ^^C + ^^C reaction 
below 150 MeV/nucleon are calculated with the antisymmetrized version of 
molecular dynamics combined with the statistical decay calculation. Density 
dependent Gogny force is used as the effective interaction. The calculated 
balance energy is about 100 MeV/nucleon, which is close to the observed 
value. Below the balance energy, the absolute value of the fragment flow is 
larger than that of nucleon flow, which is also in accordance with data. The 
dependence of the flow on the stochastic collision cross section and its origin 
are discussed. All the results are naturally understood by introducing the 
concept of two components of flow: the flow of dynamically emitted nucleons 
and the flow of the nuclear matter which contributes to both the flow of 
fragments and the flow of nucleons due to the statistical decay. 
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I. INTRODUCTION 



The collective flow in heavy ion reactions in the intermediate and high energy region 
has been studied in these years both experimentally and theoretically with the expectation 
that it carries the information of hot and dense nuclear matter. However, by the study with 
microscopic simulation approaches such as Vlasov-Uehling-Uhlenbeck (VUU) method, it has 
turned out that extracting the equation of state from the collective flow is not easy since it 
reflects not only the equation of state but also the cross section of the two-nucleon collision 
term which has theoretical ambiguity in the nuclear medium. 

Although the collective flow is a one-body observable, the most characteristic feature of 
intermediate energy heavy ion reactions is the fragment formation. Recently the collective 
flow has become to be measured with the identification of charges and/or mass numbers of 
fragments [0,0 . The large flow of fragments is observed compared to the flow of nucleons, 
which suggests that the flow of composite fragments carries the direct information of the 
equation of state (or the mean field), while most nucleons are emitted by the hard stochastic 
collisions which erase the effect of the mean field. Thus we can expect that the study of 
the fragment flow together with the nucleon flow may give us important information on the 
equation of state. 

For the theoretical analysis of the flow of fragments, the model should be able to describe 
the dynamical fragment formation, and we use the antisymmetrized version of molecular dy- 
namics (AMD) which was constructed by incorporating the stochastic collision process 
into the fermionic molecular dynamics that describes the system with a Slater determinant 
of Gaussian wave packets We have demonstrated that the AMD can describe some 

quantum mechanical features such as the shell effect in the cross sections of the dynamical 
production of the fragments. Since the flow of a particles is important among the flows 
of various fragments, this ability of the AMD to describe the large cross sections of the 
dynamical production of a particles is a very advantageous property of the AMD for the 
present study. Furthermore, AMD has an advantage that the calculation with finite-range 
effective interaction is as feasible as the calculation with zero-range effective interaction, and 
therefore the momentum dependence of the mean field, which is essential for the study of 
the equation of state, is automatically taken into the calculation. 

In studying the collective flow, we should also pay attention to the momentum distri- 
bution of the produced fragments. Namely, since the flow is known to be a very sensitive 
quantity to the various parameters in the theory, such dependence should be understood on 
the level of the momentum distribution. For example, some change of the equation of state 
and another change of the in-medium cross section may result in the same change of the 
flow value, but this does not mean that one must give up the determination of the equation 
of state by using heavy ion reactions. We still have chance that these two changes will result 
in the distinguishable changes of the momentum distribution of fragments. Our framework 
AMD is suitable for such study of the momentum distribution of fragments since there is 
no ambiguity in identifying the fragments. As we pointed out in Ref. ||^, the dissipated 
component of the momentum distribution of fragments such as a particles is sensitive to 
the in-medium cross section, while the flow angle is expected to be sensitive to the effective 
interaction. 

In this paper, we will investigate the flow of nucleons and fragments for the reaction 
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i2q_|_i2q ^Yie energy region 45 MeV/nucleon < E < 150 MeV/nucleon which includes the 
observed balance energy. After giving a brief explanation of our framework in Sect. II, we 
discuss in Sect. Ill the use of the finite range Gogny force |^ as the effective interaction which 
gives the soft equation of state with the incompressibility of the nuclear matter K = 228 
MeV. It is shown that the Gogny force reproduces binding energies of light nuclei relevant 
to the present study. In Sect. IV, before showing the calculated results of the flow, we will 
check the validity of our calculation by comparing the proton spectra with the data and the 
results of other calculations. In Sect. V, the calculated results of the flow are shown and 
it will turn out that the soft equation of state with the Gogny force is consistent to the 
experimentally observed balance energy of the flow. The calculated flow of a particles is 
larger than the flow of nucleons, which is also consistent to the feature of the experimental 
data, and the mechanism of the creation of these flows will be revealed by paying attention 
to the time scale of the reaction and dividing the collective flow into two components, i.e., 
the flow of dynamically emitted nucleons and the flow of the nuclear matter. In Sect. VI, 
the dependence of the flow on the in-medium cross section and the origin of this dependence 
are discussed. Sect. VII is devoted to the summary. 



II. FRAMEWORK 

Since the framework of the antisymmetrized version of molecular dynamics (AMD) was 
described in detail in Refs. |^,^ , here is shown only the outline of our framework. 

In AMD, the wave function of yl-nucleon system |$) is described by a Slater determinant 
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where represents the spin-isospin label of zth single particle state, Q;i = pt,pi,n|, or 
n I, and x is the spin-isospin wave function. (pZi is the spatial wave function of ith single 
particle state which is a Gaussian wave packet 



(2) 



where the width parameter u is treated as time-independent in our model. We took = 0.16 
fm~^ in the calculation presented in this paper. 

The time developments of the centers of Gaussian wave packets, {Zj(z = 1,2, . . . , A)}, 
are determined by two processes. One is the time development determined by the time- 
dependent variational principle 
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which leads to the equation of motion for {Z}. 

The second process which determines the time development of the system is the stochastic 
collision process due to the residual interaction. We incorporate this process in the similar 
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way to the quantum molecular dynamics (QMD) by introducing the physical coordinates 
{W,} li as 



where 



W. = E(\/Q)./., (4) 



Q^j = 1^77^^^^og{<!>{Z)mZ)). (5) 



d{Z* ■ Z 



In addition to the usual two nucleon collision process, the stochastic nucleon-alpha collision 
process is also included as in Ref . . Since we apply AMD to the reactions with higher inci- 
dent energy in this paper than in the previous works , the cross sections and the angular 
distributions of the stochastic collisions are reparametrized as explained in the Appendix B. 

Since in AMD the center-of-mass motion of the fragment is described by the Gaussian 
wave packet, we have to subtract from the AMD Hamiltonian (<l>(Z)|if|$(Z))/($(Z)|$(Z)) 
the sum of the spurious zero-point energies of the fragments whose number changes in 
time. The prescription to deal with this problem which is a revised version of the previous 
prescription in Refs. is explained in the Appendix C. 

The simulations of AMD are truncated at a finite time (150 fm/c in the calculation 
presented in this paper). The dynamical stage of the reaction has finished by this time and 
some excited fragments have been formed which will emit lighter particles with a long time 
scale. Such statistical decays of the equilibrated fragments are calculated with a code of 



Ref. [10 1 which is similar to the code of Piihlhofer ITT 



III. CHOICE OF THE EFFECTIVE INTERACTION 

Since the collective flow is expected to reflect the momentum dependence and the density 
dependence of the mean field, the choice of the effective interaction in our framework is 
important. 

The momentum dependence of the mean field can be taken into account as long as a finite 
range effective interaction is used, since the system is described with an antisymmetrized 
wave function in AMD. The calculation with a finite range force is as feasible as the cal- 
culation with a zero-range force in AMD, and therefore we have been using finite range 
effective interactions. In Refs. we used the Volkov force No. 1 |jl2[ with the Majorana 



parameter m = 0.576, which is a density-independent two-range interaction. By adjusting 
the parameter related to the subtraction of the spurious zero-point oscillation of fragments, 
the binding energies of nuclei lighter than ^^C are reproduced very well. However, heavier 
nuclei are overbinding with this effective interaction as shown in Fig. |l|. When the Volkov 
force is applied, usually its Majorana parameter m should be chosen depending on the mass 
number of the nuclei under consideration, which makes the Volkov force very inconvenient 
for our study. 

In order to reproduce the binding energies of nuclei in wide mass number region within a 
single effective interaction, and in order to cope duely with the formation of the high density 
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FIG. 1. Binding energies per nucleon of nuclei calculated with the Volkov force No. 1 
(m = 0.576) (left) and the Gogny force (right). Isotopes are connected with lines. Experimental 
data shifted by —5 MeV are also shown. 



nuclear matter in collisions with high incident energy around and above 100 MeV/nucleon, 
the density dependence of the effective interaction should be introduced. We use the Gogny 
force 1^ in the calculation presented in this paper. The binding energies of nuclei lighter 
than 2^Mg are reproduced (Fig. |l|) by adjusting the parameter related to the subtraction 
of the spurious zero-point oscillation of fragments which is described in the Appendix C. It 
is known that the Gogny force can describe the saturation property of nuclear matter with 
the incompressibility K = 228 MeV. The momentum dependence of the mean field with the 
Gogny force is reliable below 200 MeV, since the energy dependence of the real part of the 
nucleon optical potential is reproduced in this energy range. The Gogny force has ever been 



successfully applied to the problem of the flow in the framework of VUU [|l^] . 

The calculation with density dependent zero-range force requires us to evaluate the quan- 
tities such as 



/ drp{r] 



2+<7 



(6) 



where p(r) is the density calculated with the AMD wave function. This spatial integration 
was performed by a kind of Monte Carlo method with test particles generated randomly 
with a weight function. The detail of the calculation method is explained in the Appendix 
A. If we use lOOA test particles, the statistical error in the evaluation of the energy of ^^C 
ground state is about 2 MeV, which seems sufficiently small for the current purpose. 
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FIG. 2. Energy spectra of protons for the reaction ^^C + ^^C at 84 MeV/nucleon. Dashed 
histograms are the proton spectra calculated before the statistical decay, and the solid histograms 
are the final calculated results. Experimental data are shown by crosses. Energy E and angles are 
in the laboratory system. 



IV. ENERGY SPECTRA OF PROTONS 



data 



Before going into the problem of the collective flow, we compare our results on the proton 
energy spectrum to the data in the reaction ^^C + ^^C reaction at the incident energy 84 
MeV/nucleon, which was also analyzed by other simulation methods such as VUU and QMD 
||14|| . It is for the purpose of showing the reliability of the AMD calculation in the incident 
energy region around 100 MeV/nucleon. As mentioned before the AMD has proved to 
be very successful in reproducing the nuclide distribution and the fragment energy spectra 
in the incident energy region around the Fermi energy. 

In Fig. 1^, the calculated proton energy spectra for various angles are compared to the 
TB|. The dotted histograms are the calculated proton spectra at the end of AMD 
simulation t = 150 fm/c, which do not contain the contribution from the statistical decay of 
excited fragments after t = 150 fm/c, while the solid histograms are the final results which 
should be compared to the data. We can see the effect of the statistical decay of fragments 
on these spectra is not so large except for the low energy part. Although our calculation 
has tendency to overestimate the proton spectra at 35°-65°, this feature and the degree of 
the overestimation are common with other simulation approaches such as VUU and QMD 
I^ . The spectra of backward angles are reproduced by our model as well as by QMD. In 
Fig. 1^, the multiplicity of nucleons are shown as a function of the impact parameter. The 
dynamical contribution before the statistical decay is shown as well as the final multiplicity. 



This result is again quite similar to those of QMD and VUU for all the impact parameter 
region. 

In calculating the spectra in our model, we must attribute some widths to the momenta of 
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FIG. 3. Impact parameter dependence of the nucleon multiplicity for the reaction ^^C + ^^C 
at 84 MeV/nucleon. Dashed histogram shows the multiplicity before the statistical decay and the 
solid histogram shows the final calculated result. 



nucleons and fragments produced in the AMD calculation, since these momenta are central 
values of Gaussian wave packets. Each wave packet has the width Ap/ i/Ap in the momentum 
per nucleon where Ap = and Ay is the mass number of the fragment. However, we 
should not always take this width to have physical meaning, since some part of this width 
comes from the unphysical width of the projectile or the target in the initial state of the 
simulation. It seems reasonable to consider that such unphysical contribution mainly goes 
to the momentum widths of fragments and to assume that all the widths of the wave packets 
of nucleons are physical. When the proton spectra in Fig. Q are calculated, we neglect the 
widths of the fragments produced dynamically in the AMD simulation, but attribute the 
widths to the dynamically produced nucleons in the AMD simulation. 



V. CALCULATED RESULTS OF THE FLOW AND THE MECHANISM OF THE 

FLOW CREATION 

Although the definition of the flow is somewhat complicated in the experimental situation 
where the determination of the reaction plane is not trivial, we simply define here the flow 
of fragments with mass number as the mean in-plane transverse momentum, 

(^^'.■)Mf, (7) 

where w = 1 if > and w = — 1 if < 0. Pj, and P^ are the components of the in- 
plane momentum of the fragment in the center-of-mass system which are perpendicular to 
and along with the beam direction respectively. Negative flow corresponds to the attractive 
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FIG. 4. Calculated flows of nucleons and a particles for various incident energies for the 



reaction ^^C + ^'^C. The impact parameter of the simulation is fixed to be 2 fm. 
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interaction between the projectile and the target. This definition of the fiow was already 
used in Ref. [0] and has almost no uncertainty due to the interpretation of the momentum 
width of the wave packets discussed in the previous section. 

Another definition of the fiow as the slope of {Px)-Pz curve is of course possible, but 
we have checked that the qualitative features to be discussed below does not change. The 
statistical precision of the fiow defined by Eq. is better than the latter one. 

Fig. H shows the energy dependence of the calculated fiows of nucleons and a particles. 
In this figure, as well as in all the results hereafter, the impact parameter of the reaction is 
fixed to be 2 fm in order to improve the statistical precision. We can see that the calculated 
fiows of nucleons and a particles change their sign at the balance energy i^bai = (100 ± 20) 
MeV, though the balance energy of a particles seems to be larger than that of nucleons. This 
calculated balance energy is not different so much from the experimentally deduced value 
122 ± 12 MeV/nucleon [^. On the other hand, in the calculation with the Volkov force, the 
fiows of nucleons and a particles are negative even at the incident energy 150 MeV/nucleon 
||T7| . (In this calculation, another version of the stochastic collision process adopted in Ref. 
P was used.) 

A remarkable feature is that the absolute value of the fiow of a particles is much larger 
than that of nucleons in the energy region E < -Ebai- The reason of this larger fiow of 
fragments becomes clear by distinguishing nucleons and fragments according to the time they 
are produced. From Fig. ^ which shows the distribution of the emission time of nucleons, we 
can see that the dynamical stage of the reaction has already finished at t = 150 fm/c. Fig. |^ 
shows the fiows att = 150 fm/c (Dynamical) and the fiows of nucleons and a particles which 
are produced by the statistical decay of other fragments after t = 150 fm/c (Evaporation). 
Although the 'dynamical' fiows in the figure contain some contribution from the equilibrated 
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FIG. 5. The dynamical flows of nucleons and a particles which are calculated before the 
statistical decay, and the flows of evaporated particles which are calculated only from the decay 
products of the statistical decay. 



stage of the reaction, we can see clear deference between the flow of dynamically emitted 
nucleons and the flow of evaporated nucleons from the equilibrated fragments. The negative 
flow of the evaporated nucleons is much larger than the flow of the dynamically emitted 
nucleons, and the final result is between them. Since larger number of nucleons are emitted 
in the dynamical stage, as can be seen from Fig. |^, the final result is closer to the flow of 
the dynamically emitted nucleons. On the other hand, the flow of a particles seems to be 
almost independent of the time they are produced, and the value of the flow is close to the 
flow of the evaporated nucleons. 

These results lead us almost uniquely to the following interpretation. We can separate 
the collective flow into two components. The first one is the flow of the nucleons which 
are emitted by the hard stochastic collisions in the dynamical stage of the reaction. Since 
the stochastic collision erases the memory of the mean field which the nucleon has been 
feeling before it is scattered, the absolute value of the negative flow is small. There may be 
another geometrical reason for this small flow that the projectile and the target prevent the 
nucleons with negative flow from going out to the free space more than the nucleons with 
positive flow. The second component of the collective flow is the flow of the nuclear matter 
which contributes to the flows of composite fragments produced in the dynamical stage of 
the reaction. This component has been affected directly by the mean field and the absolute 
value of the negative flow is large. The flows of evaporated nucleons and fragments inherit 
the flow of the parent fragments and hence they have the large value similar to the original 
nuclear matter flow. 
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FIG. 6. Time dependence of the nucleon emission rate for the central reaction ^^C + ^^C 
(5 = 2 fm) at various incident energies. A nucleon is interpreted to have been emitted when no 
other nucleons exist within a radius of 3 fm. 
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FIG. 8. Dependence of the flows on the stochastic colhsion cross sections for the reaction 
12q _^ 12q 45 MeV/nucleon. 

VI. DEPENDENCE ON THE STOCHASTIC COLLISION CROSS SECTION 



Since the stochastic colhsion cross sections have some theoretical uncertainty, we should 
investigate the dependence of the flow on them. Furthermore the study of such a-dependence 
of the momentum distribution will give us the way to determine the in-medium cross sections. 

Figs. H and |^ show the cx-dependence of the flows of nucleons and a particles at the 
incident energy 45 MeV/nucleon and 110 MeV/nucleon, respectively, when all the cross 
sections (see Appendix B) are multiplied by 0.5, 1.0 and 1.5. The flows of the dynamically 
produced particles and the flows of the particles produced by the statistical decay are shown 
together with the flnal results. The flow of the dynamically emitted nucleons is almost 
independent of the cross sections. The multiplicity of dynamically emitted nucleons in the 
reaction at 45 MeV/nucleon is 2.00, 3.52 and 4.95 for 0.5cr, l.Ocr and 1.5cr, respectively. 
These results are consistent with the interpretation in the previous section because what is 
important to the flow of the dynamically emitted nucleons is that the emitted nucleons are 
direct products of the stochastic collisions and the number of such nucleons is not relevant 
to the value of the flow. 

On the other hand, the cx-dependence of the flow of the nuclear matter is large, especially 
in the reaction at 45 MeV/nucleon, and the absolute value of the negative flow is smaller 
for larger cross sections. This might sound inconsistent to the naive expectation from the 
previous interpretation that the flow of the matter is essentially affected by the mean fleld. 
This cr- dependence, however, turns out to be very reasonable when one pays attention to 



other features of the momentum distribution as well as the flow. Figs. |T0| and [11] show the 
in-plane momentum distribution of nucleons and a particles which are calculated at t = 150 
fm/c (Dynamical) and from those produced by the statistical decay process (Evaporation), 



11 



^^C+^^C b=2fm 



> 





< 

A 
Q. 



30 



20 



10 



-10 



-20 



-30 



1 1 MeV/u 
Nucleon 



Final 
Dynamical 
Evaporation 




> 





< 

X 

Q. 

V 



30 



20 



10 



-10 



-20 



-30 



1 1 MeV/u 
Alpha 


Final <> 
Dynamical 
Evaporation q 




Gogny 


0.5 


1 1.5 


Cross Section [a] 



0.5 1 1.5 

Cross Section [a] 

FIG. 9. Dependence of the flows on the stochastic colhsion cross sections for the reaction 
12q _^ 12q ;^;^Q MeV/nucleon. 



respectively. Results for three cross sections O.So", l.Oo" and l.So" are shown for the reaction 
at 45 MeV/nucleon. The momentum distribution of dynamically emitted nucleons is almost 
spherical while the other three momentum distributions (of evaporated nucleons, dynamical 
a particles, and evaporated a particles) have the common cr-dependence though the mo- 
mentum distribution of the evaporated nucleons has wider spreading compared to that of a 
particles. As the cross sections are increased, the dissipated component of the momentum 
distribution around the center-of-mass momentum increases while the flow angle does not 
change. What has changed the flow value of the matter is nothing but this cr-dependence of 
the dissipated component of the momentum distribution. 

It is to be noted here that the cx-dependence of the flow is not due to the specific definition 
of the flow (0). For example, if we adopt the definition of the flow as the slope of {Px)-Pz 
curve at P2 = 0, we get the cr-dependence of the nucleon flow shown in Fig. |1^ which has 
the same qualitative feature as the left part of Fig. |. 



VII. SUMMARY 

In this paper the collective flow of fragments was calculated as well as the flow of nucleons 
with the AMD for the reaction ^^C + ^^C in the energy range 45 MeV/nucleon < E < 150 
MeV/nucleon. The balance energy observed in the experiment was found to be consistent 
with the Gogny force which yields a soft equation of state with K = 228 MeV. We obtained 
larger absolute value of the fragment flow than the nucleon flow, which is also consistent 
with the experimentally observed feature. 

The energy dependence and the cr-dependence of the calculated flows of nucleons and 
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FIG. 10. Momentum distributions of nucleons and a particles at t = 150 fm/c before the 
statistical decay in the reaction ^^C + ^^C at 45 MeV/nucleon with 6 = 2 fm. Results for the three 
sets 0.5(7, 1.0(7 and 1.5a of the stochastic collision cross sections are shown. 
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fragments were found to be understandable, at least in the energy region where the flow is 
negative, by introducing the idea of two components of the collective flow. The first com- 
ponent is the flow of the dynamically emitted nucleons which are produced by stochastic 
collisions in the early stage of the reaction. Since the stochastic collisions erase the effect of 
the mean field, the fiow of the dynamically emitted nucleons is small. The second compo- 
nent is the fiow of the nuclear matter which contributes to the fiows of fragments produced 
during the dynamical stage of the reaction. The flows of the evaporated nucleons and frag- 
ments inherit the flow of the parent fragments and hence the flows of dynamical fragments 
and the evaporated nucleons and fragments have almost common magnitude and common 
energy dependence and o"-dependence. Although the second component is affected largely 
by the mean field, its cr-dependence is also large. However, this cr-dependence turned out to 
come not from the change of the flow angle but from change of the yield of the dissipated 
component in the momentum distribution of the fragments. Detailed comparison of the 
fragment momentum distribution with the data will be able to settle the in-medium cross 
sections. On the other hand, nucleon momentum distribution is mainly governed by the dy- 
namically emitted nucleons in this energy region, and therefore the shape of the momentum 
distribution is not so sensitive to the in-medium cross sections. 

We have mainly concentrated on the mechanism of the flow creation and the cr- 
dependence in this paper, but the dependence on the effective interaction was not discussed 
in detail. What we have found is that the Gogny force is sufficient for the reproduction of 
the balance energy while the Volkov force fails. To calculate the fiow with other effective 
interactions which reproduce the saturation property of the nuclear matter and have various 
incompressibilities is an important future problem for the determination of the equation of 
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state of the nuclear matter. 
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APPENDIX A: DENSITY DEPENDENT INTERACTION 

Here we consider an effective interaction which contains density dependent zero-range 
part, 

v, = ^-^Y.x^,P''imn~f,). (Ai) 

i<j 

For example, the density dependent part of the Gogny force has this form with 

= 6t3, X = l + P^, a = 1/3. (A2) 

By using the density p{r) which is calculated from the AMD wave function, the expectation 
value of Vp can be written as 

Vp=(V;) = ^/rfr>(r>'^(r), (A3) 



where 



/^(^ = IIPa(^Pc.(r), p{r) = Y,pa{r), (A4) 

a a 



Pa\' j \ \ ^ ^ ^13-^11 ' 

Po.{^ = E - /5«)P/3(r). (A6) 

Since the analytical integration of ( |A3| ) is impossible for non-integer a and even in the case 
of cr = 1 the evaluation of the analytical six fold summation is not feasible, we evaluate Vp 
by the Monte Carlo method generating Atp test particles randomly with an appropriate 
weight function Pw{r), 

^ E ^P'^lr.), F= [ dfp„{r), (A7) 

16 A^TP ^1 Pw(^fc) J 

where is the position of the kth test particle. The arbitrary function p^{f) should be 
chosen so as to cover the important region in the integral ([A3|) efficiently and at the same 
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time to make the generation of test particles numerically easy. By using the approximate 
density Pw(0 calculated from the physical positions Ri = ReWi as 

pArl = {-f'j:e~''^'-'^'\ (A8) 

we take /iw(^ as 



^^^i^y/^^^.0^.n^f_ (AlO) 

Derivatives of Vp with respect to are necessary to solve the equation of motion of 
AMD. They are also evaluated by the Monte Carlo method as 



dVp tp 



dr 



2p«(0 + ^1p'(.^%£^ (All) 



dZl 16 y L p(f 

leA/^TP^i^ p(rfc) Vw(rfc) dzi 

where a is the spin-isospin label of the nucleon h. 

APPENDIX B: DETAIL OF THE STOCHASTIC COLLISION PROCESS 

In the higher energy region, we use the proton-neutron and proton-proton collision cross 
sections parametrized as 



rp„ = max{l3335E[MeV]"^■^^^40}mb, (Bl) 
= max{4445 E[MeV]~^•^^^ 25} mb, (B2) 



where E is the energy in the laboratory system of the stochastic collisions. These are based 
on the data of free cross sections. The neutron-neutron collision cross section is assumed to 
be identical to the proton-proton collision cross section. The nucleon-alpha collision total 
cross section aNa,tot is chosen to be equal to the total cross section of the process where the 
nucleon were scattered by two-nucleon collisions with four nucleons in the a cluster with 
the cross sections cXp^ and a pp. The adopted nucleon-alpha inelastic collision cross section is 
based on the experimental data parametrized as 

(TNaMeX = maX {l20 - 162e-(^-2°^<=^)/(i°^<=^\0}mb. (B3) 
The angular distribution of proton- neutron scatterings are taken as 

^ oc io-"(-/2-|e-V2|)^ ^ = -max|0.3331ni?[MeV] - 1, o|, (B4) 
ail 71 ^ J 
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while the proton-proton and neutron-neutron scatterings are assumed to be isotropic. The 
nucleon-alpha inelastic collision is treated as a two-nucleon collision between the nucleon and 
a nucleon in the a cluster and the same angular distribution as the two-nucleon collision is 
applied. The angular distribution of the nucleon-alpha elastic scattering is parametrized @ 

as 



o,el 



oc exp 



180° - d^i 



70° 



+ 10 exp 



d - 20°x2 



40° 



(B5) 



Since the above defined cross sections are too large if they are applied to low energy 
region, we use the cross sections 



100 mb 



O". 



VV 



0'Na,tot 



1 + E/(200MeV) + 2min{(p/po)^/^ 1}' 
571 mb 

1 + E/(200 MeV) + 2 min{(p/po)^/^ 1} 



(B6) 
(B7) 



which are similar to those we used in Ref. for the reaction at 28.7 MeV/nucleon. The 
cross sections ( P6|) and ( P7|) are used if they are smaller than the cross sections (Tp„, app and 
c"Afa,tot defined for the higher energy region by Eqs. (|B1D and (|B^ ). On the other hand, the 
same expressions of angular distributions and cr Ara,inei are used for the whole energy region. 



APPENDIX C: SUBTRACTION OF SPURIOUS ZERO-POINT OSCILLATION 



Like in Refs. 0-|5| the spurious kinetic energies of zero-point oscillations of center-of-mass 
motion of fragments are subtracted by modifying the Hamiltonian in the equation of motion 
of AMD as 

n={H)-^A + n{A~N^), (CI) 

where the zero-point kinetic energy Tq is treated as a free parameter to adjust the gross 
feature of the binding energies of nuclei. 
We define the fragment number Np as 

iVp = E (C2) 



where 



A A 



Hi 



E — E ~f^j-> ~ E /'i (^^) 



and 



fij = F{dij,i,d), fij = F{dij,^,a), fij = F{dij,^,d), (C4) 
d^, = \ReiZi - Z,)l (C5) 
^, , . , f 1 if d < a 



a 
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TABLE I. Parameters concerned with the subtraction of spurious zero-point oscillations of 
fragments. Parameters which are used with the Gogny force in this paper are shown as well as 
those used with the Volkov force No. 1 with m = 0.576. 



Force 




a 


i 


a 




a 


go 


a 


M 


To 


Volkov 


1.0 


0.2 


2.0 


0.1 












7.7 MeV 


Gogny 


2.0 


0.6 


2.0 


0.2 


1.0 


0.5 


1.0 


2.0 


12.0 


9.2 MeV 



Compared to Ref. 0] where the Volkov force was used as the effective interaction, we have 
newly introduced a function 

g{k) = l+goe~^''''^'/"^\ (C7) 

which should be equal to 1 in principle. We have introduced it in order to remedy the 
situation that nuclei around ^^C are underbound while the binding energies of other nuclei 
such as a and ^^O are almost reproduced in the calculation with the Gogny force. We have 
chosen the parameters as shown in Table |. Binding energies of nuclei are sensitive only to 
two parameters Tq and M. 
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